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ABSTRACT 

In the context of Computer Aided Diagnosis system for di¬ 
abetic retinopathy, we present a novel method for detection 
of exudates and their classification for disease severity pre¬ 
diction. The method is based on Gaussian scale space based 
interest map and mathematical morphology. lit makes use 
of support vector machine for classification and location 
information of the optic disc and the macula region for sever¬ 
ity prediction. It can efficiently handle luminance variation 
and it is suitable for varied sized exudates. The method 
has been probed in publicly available DIARETDB1V2 and 
e-ophthaEX databases. Eor exudate detection the proposed 
method achieved a sensitivity of 96.54% and prediction of 
98.35% in DIARETDB1V2 database. 

Index Terms — Exudate, Diabetic Rationopathy, Image 
Processing 

1. INTRODUCTION 

In recent days diabetic retinopathy (DR) is one of the se¬ 
vere eye diseases causing blindness. With early stage de¬ 
tection and treatment the patient can be saved from losing 
sight. Automatic computer aided diagnosis system will re¬ 
duce burden on specialists. Also for monitoring and check¬ 
ing the progress of disease efficiently, automatic system will 
perform much better than human in terms of manual evalua¬ 
tion time. Since comparison and evaluation of images manu¬ 
ally is a time consuming task and images are subject to var¬ 
ious distortions. Eor accurate analysis of progress of dia¬ 
betic retinopathy, detection of exudate is mandatory. Exu¬ 
dates are primary clinical symptoms of diabetic retinopathy. 
Two types of exudates namely soft exudate and hard exudate 
appear. Hard exudates are visible in non-proliferative diabetic 
retinopathy and soft exudates (cotton wool spots) in prolifer¬ 
ative diabetic retinopathy. Hard exudates represent the accu¬ 
mulation of lipid in or under the retina secondary to vascular 
leakage and visible as discrete yellowish deposits in color fun¬ 
dus images. Cotton-wool spots are nerve fibre layer infarcts 
and they are visible as pale white rather than yellow. Also 
exudates are variable in sizes and shapes. A typical patholog¬ 
ical retinal image is depicted in Eig. 1 to show features like 
the optic disc, the macula, the blood vessels and exudates. 
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Fig. 1: Retinal Eeatures 


Eor identifying the stage of DR, classification of soft and 
hard exudates is foremost important to distinguish them for 
other retinal pathological features like drusen, heamorrhages, 
microaneurysms etc. Eundus images are prone to artifacts 
related to defocus, motion blur, fingerprints etc. This arti¬ 
facts are prone to interpretation as pathological features be¬ 
cause of their similarity with exudates and drusen. One im¬ 
portant defect of fundus image is luminosity and contrast vari¬ 
ation, which is improperly addressed in many exudate detec¬ 
tion methods. 

Several methods have been presented for detection of ex¬ 
udates in colour fundus photograph using image processing 
and machine learning algorithms. Sanchez et al. O use a 
mixture model to separate exudate from background. An edge 
detection based method was used to remove other outliers. 
Giancardo et al. Q use a Kirsch’s edge method to assign 
score for exudate candidate to a pre-processed . They have 
used background estimation and image normalization for pre¬ 
processing. Eor classification and detection of drusen, exu¬ 
dates and cotton wool spots, a pixel wise classification algo¬ 
rithm is presented by Niemeijer et al. a. Zhang et al. SI use 
mathematical morphological and contextual features for can¬ 
didate extraction followed by random forest based classifica¬ 
tion method for exudate detection. Euzzy c-means clustering 




based method was used for segmentation of features and neu¬ 
ral network classifier for exudate detection by Osareh et al. 
0. Rocha etal. El have addressed the problem of detect¬ 
ing bright and red lession by using one novel algorithm, they 
made use of SURF features with machine learning method. 
But they failed to properly address the problems of luminance 
variation and artifacts. 

In this work,we propose a novel exudate detection method 
using Gaussian scale space based interest map (GIMAP) and 
mathematical morphology. This approach is robust to arte¬ 
facts and illuminance variation. Secondly a disease severity 
prediction method is developed by using information of exu¬ 
date location with respect to the macula region and the optic 
disc. In addition to that we propose classification system us¬ 
ing SVM for hard and soft exudate. Distinguishing exudate as 
hard and soft important for severity prediction also to identify 
type of diabetic retinopathy, whether it is non-proliferative or 
proliferative. Section 2 presents our method. Experimental 
setup and the results are discussed in section 3 and finally 
conclusions are are drawn in section 4. 

2. METHOD 
2.1. Exudate Detection 

Exudate are the bright lessons found in retinal image, caused 
due to diabetic retinopathy, a most common disorder of eye 
with patient having diabetes. It is also a main reason of blind¬ 
ness. For detection of the optic disc and fovea, the method 
described by Niemeijer et al. ||6| was used. Steps involve in 
this method is shown in Fig. 2. 
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Fig. 2 : Method Overview 

which algorithm will be iterated to find solution. The diffu¬ 
sion function c(r, 6, a) is a monotonically decreasing func¬ 
tion of the image gradient magnitude. 


2.7.7. Preprocessing 

For correct detection and classification of hard and soft 
exudates it is important to reduce noise while preserving 
edges. In this method anisotropic diffusion filtering |l7l was 
used for reducing noise in the images and preserving edges. 
Anisotropic diffusion filters have been proved to be succesful 
in edge preserving smooting and denoising of medical im¬ 
ages. Hard and soft exudates possesses discriminative edge 
structures, preserving this edge while reducing noise is also 
one of important step towards classification. Step involves in 
filtering is as follows (Tl. 


■^I{r,g,b,a) = V • (c{r,g,b,a)\/I{r,g,b,a)) (1) 


c{r,g,b,a) = f{\VI{r,g,b,a)\) (2) 

= 1 I > 0 

where I is the input image, (r,g,b) is color channel of it, K 
is diffusion constant and a is noise standard deviation over 


2.7.2. Gaussian Scale Space Construction 

Gaussian Scale space based interest map can capture local 
image structure and scale invariant image features. Since ex¬ 
udate varies on the basis of size and hence different exudate 
will response to different scales of Gaussian. For identify¬ 
ing all exudate interest map using Gaussian scale space was 
constructed M- First step in the construction of GIMAP is 
the computation of 1st derivative of Gaussian as filters at sev¬ 
eral scales and smoothed the derivative using Gaussian filter 
and take absolute values of derivatives. In addition to that 
Laplacian of Gaussian for each scale was computed. Because 
of colour difference between hard and soft exudates, they re¬ 
spond variably to colour channels. Now taking the maximum 
response of all the colour channels final interest map for the 
scale is constructed. This process is repeated over different 
scales. For each scale we have two filter 1st derivative Gaus¬ 
sian and Laplacian of Gaussian, take maximum of absolute 
values from both filter ouputs in the interest map. Since our 
concerned features may appear in variable scales so interest 
map of each individual map is combined using maximum op¬ 
eration to form a decision making interest map. Scale used in 
these works are (\/2, h/2...k\/2). The decision making inter- 













est map also contain some outliers such blood vessels, vein, 
haemorrhage etc. Filtering process is described by following 
equations ( 91 . 


^nd cr{x,y) be the mean and standard deviation of win¬ 
dow cantered at (x, y), then threshold th{x, y) will be defined 
as follows. 
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For computation of Laplacian of gaussian we have used dif¬ 
ference of Gaussian approximation. If I{v, k, f) be a filtered 
image, v[r^g^ b] is its three color channel at a scale k using 
filter /, then the interest map for this scale is obtained by se¬ 
lecting maximum response over color channel as follows. 


I{k,f) = maxyl {v,k,f) (7) 


And the decision making interest map is obtained taking max¬ 
imum over filter / and then scale k as follows. 


=TO(^,y)[l + c(^^^ - 1)] (12) 

where S is maximum of standard deviation of all windows 
and c G [0.20.5] is a parameter. 

Along with exudate other features such as blood vessels 
and veins is also detected due to Gaussian scale space. Mor¬ 
phological opening operation and connected component anal¬ 
ysis was used for removing those unwanted features. Blood 
vessels and veins are characterized by thin long geometrical 
structures. A rectangle structural element having width and 
length larger than the thickness of those blood vessels is used 
for removing those structures. Remaining vessels and outliers 
left after this operation will be removed by using concept of 
convexity and connected component analysis. For each con¬ 
nected region we will compute convex hull of the region and 
define the following prohReg eq. (13) term. 


prohReg = 


^i,jGH 

YjijR 


(13) 


dMap = maxkmaxfl(k^ f) (8) 

2.1.3. Binarization and Vessel Removal 

The decision map obtained from above procedure need to be 
improved by reducing outliers present. First step is to en¬ 
hance the features by using grayscale morphological opera¬ 
tions. Here closing operation with a ’disk’ structuring element 
of size 2 and 3 was used. 

Dilation : {f ^ B){x^y) = maxf{x — s^y — t)\{sR) G B 

Erosion : {f Q B){x,y) = maxf{x + 5, ^ + t)|(s, t) G B 

(9) 

Opening : {f o B) = {f e B) ® B 

Closing : {f • B) = (f ® B) e B ^ 


The prohReg term computes probability of a connected re¬ 
gion being exudate or not. It is the ratio of the area of the re¬ 
gion convex hull and its area. This term will give approximate 
idea of the connected region. If the value of prohReg >0.8 
then the region will be discarded, the connected region will 
be considered as vessel structures or fiash artifacts. 

where GH is convex hull of the region R and (i, j) de¬ 
notes the pixel positions. If convexity of a connected region 
is less than 0.2 then the area will be discarded. Also we want 
to ensure that artifacts like fiares gets removed from GIMAR 
From observation it has been noticed that fiares are circular 
shape object. None of the features like hard and soft exudate 
have shapes like that of fiares, approximate circular. Com¬ 
pactness (G) is defined below eq. (14) of every connected 
region of the GIMAP should be less than J. 


If iMapl and iMap2 are obtained after closing opera¬ 
tion, then resulting enhanced image is obatined by following 
operation. 


n 


(14) 


iMap = maxij {iMapl GMap2) (11) 

where (i, j) are pixels position in both images. 

In the second step convert the interest map to binary map. 
For binarization of the interest map Sauvola’s Col local adap¬ 
tive thresholding technique is used. Local threshholding is 
efficient in this particular situation beacuse pixels values of 
hard and soft exudates vary significantly and to get both of 
those in final map we need to binarize using local windows. 
Let /(x, y) is interest map image, take a window size of 9 x 9, 


where Ar and are area and perimeter of the region R re¬ 
spectively. In the final step we will use colour and luminance 
features to select exudates and remove outliers, in general we 
didn’t observe outliers like fiares, flash etc after adapting pre- 
vios mentioned steps. But still there are possibility fo includ¬ 
ing drusens. Hard and soft exudates are characterized by its 
yellowish and whitish color significantly different from that 
of haemorhages, microanerysms and other related features. 
At the end, a support vector machine classifier is used to get 
final labels of hard and soft exudates. 












2.2. Classification of Exudates using S VM 

Two types of exudates namely hard and soft are associated 
with different diabetic retinopathy symptoms. We will use 
a Support Vector Machine classifier to differentiate between 
soft and hard exudates and outliers. SVM is one of most effec¬ 
tive classifier, with its good generalization ability. Selection 
and extraction of feature vector is one of the challenging task 
for efficient classification. Since hard and soft exudates are 
varied in color magnitude we will use mean R, G, B values 
as feature for each detected regions. In addition, to ensure 
luminance invariance in classification process, images will be 
converted to Lab color space and values of a, b channels will 
be used as features. Both hard and soft exudates are varied in 
terms of shape and edge structures. Hard exudates are chrac- 
terizes by uneven edge and soft exudates are by smooth edge 
and circular structures. Firstly we will design features vector 
for classification by using following data. Label 1 shows the 
features used in this method for SVM. 


Table 1: Features for SVM 


Feature 

Description 

1 

Area of Connected Region 

2-4 

Mean R, G, B Intensity Values 

5-7 

Standard deviation of R, G, B intensity values 

8-9 

Mean a, b intensty values in Lab color space 

10-11 

Standard deviation of a,b intensity values 

12 

Ecentricity of the region 

13 

Extent of the region 

14-15 

Major and Minor axis length 

16 

Convexity 

17 

Gradient values of edge pixels of the region 

18 

Compactness of the Region 

19 

Energy of the Region 

20-22 

Color contrast with neighbour regions 


Features vector constructed from above ideas is 22 di¬ 
mensional. For SVM classification, we have three class 
specifically hard exudate, soft exudate and outliers. System is 
trained with hard, soft exudates and non exudates pixels data 
including artifacts from expert labeled images. 

2.3. Severity Detection 

As per ophthalmologists, exudate (or other retinal pathology 
such as microaneurysms, heamorrhages, drusen etc.) location 
with respect to the fovea and the optic disc determine severity 
of the diseases. If exudates location are close to those features 
and causes defect to those features then patients may lose vi¬ 
sion. As per the norm of the International Council of Oph¬ 
thalmology Diabetic Retinopathy IT4l divided in five groups 
specifically: none, mild, moderate, severe, and proliferate. 
We will define some rule to detect severity by using the fovea 



Fig. 3: Severity level Groups 


and the optic disc location as centre of two set of concentric 
circles. Radius of each concentric circles with the fovea lo¬ 
cation as centre is multiple of 80 pixels. The radius selection 
for each circle based on the average value of the optic disc 
width and height. Also for each concentric circles with the 
optic disc centre location as circle centre is multiple of 55 pix¬ 
els. This setting valid for images with width 1500 pixels and 
height 1152 pixels. If exudates reside in vey inner circle then 
we label it as Leveli, other consecutive circle will be known 
as Levels, where N is the index number of the circle, bigger 
N means circle is of bigger radius. For each circle area of ex- 
udate(No of pixels) will be calculated. Finally by using area 
and location information diseases severity will be predicted. 
If location of exudate in Leveli circle then on the basis of 
its areas it will included in proliferate or severe groups. Here 
area of Levels circle means area of {Levels — Leveljsf-i) 
circle. Fig. 3 depicts the scenarios described above. 

Severity oc (ci * A + C 2 * —) 

c. = m ( 15 ) 

C2 = f{D) 

where A denotes area of exudate region, D its distance from 
either the optic disc center or the fovea and ci, C 2 G [0,1] are 
two constant function of distance. 

Proliferate : Xi G Ci D > -—Ac^ (16) 

16 

Severe : (xi € Ci fl Y,iXi < 

(17) 

{xi e C2 n Y,iXi > y^^c2) 

Moderate : {xi e (72 fl Sja;, < —Ac^)\J 

(18) 

{xi e C2^T.iXi > —Ac^) 
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Fig. 4: Left: Original Image, Right: GIMAP, Down - Left: 
After Binarization, Down - Right: Final detected Exudate 


Mild : {xi G C3 n T.iXi < —Ac3)U 


{xi G C2 nSjXj > 


(19) 


where Xi is a exudate pixel, Ci denotes Leveli circle and 
Aci its corresponding area. 


3. RESULT AND DISCUSSIONS 

For analysing the accuracy of our method on different set 
of images, we have used publicly available DIARETDBl 
CD and e-ophthaEX oa dataset. Images with various diffi¬ 
culty from these databases have been chosen for testing our 
algorithm accuracy, this dataset include image with patho¬ 
logical features such as haemorrhages, exudates and microa¬ 
neurysms. This dataset provided by experts in Ophthalmol¬ 
ogy with proper pixel wise annotation of features location 
in images. This three dataset includes 400 retinal images 
with variety of pathological symptoms. Eor implemetation 
MATLAB platform in a windows 8.1 machine with Intel i7 
processor was used. 

Eor detection of exudate Gaussian scale space was con¬ 
structed using 10 different scale. At first each images is re¬ 
sized by a factor of- ? - —j -^ using cubic inter- 

polation method for reducing computational time. Eor each 
scale first derivative of Gaussian and Laplacian of Gaussian is 
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Fig. 5: Sensitivity and prediction variation 


computed for analysing structure of exudates present in reti¬ 
nal images. By using those values a decision making interest 
map is formed. Different types of outliers were removed us¬ 
ing morphological connected component analysis and open¬ 
ing operation. 

Eor accuracy analysis of exudates detection we will com¬ 
pute true positive (TP) a number of exudates pixels correctly 
detected, false positive (EP) a number of non-exudate pixels 
which are detected wrongly as exudate pixels, false negative 
(EN) number of exudate pixels that were not detected and true 
negative (TN) a number of no exudates pixels which were cor¬ 
rectly identified as non-exudate pixels. Also sensitivity and 
specificity at pixel level is computed. Thus the global sensiv- 
ity SE and the global specificty SP and accuracy AC for each 
image is defined as follows. 


AC = 


SE = 


TP 


tptfn 

TP 

tptfp 

TN 

TNPFP 
TP PTN 

TPPTNPFPT FN 


PRFD = 


SP = 


( 20 ) 


A detailed result of accuracy obtained is illustrated in Ta¬ 
ble 1. 


Table 2: Result of Exudate detection 


Database 

Resolution 

Sensitivity Prediction AUC 

eophthaEX 

960 X 1440 

95.82 

97.85 

0.962 

DIARETDBlv21152 x 1500 

96.54 

98.35 

0.968 


Eig. 4 shows exudate detection intermediate and final step 
results. Eirst row shows the original image and its correspond¬ 
ing GIMAP, where as in second row image obtained after bi¬ 
narization and final detection result is depicted. Prom this 
Pigure it can be noticed that GIMAP respones to retinal fea- 










































Table 3: Result of Exudate Classification 


Database 

Resolution 

# images 

Type 

True PositivePer image 

Ealse PositivePer image 

eophthaEX 

960 X 1440 

80 

Hard Exudate 

98.73 % 

3.05 % 

DIARETDBlv2 

1152 X 1500 

80 

Hard Exudate 

98.51 % 

2.16% 

DIARETDBlv2 

1152 X 1500 

80 

Soft Exudate 

98.23 % 

2.35 % 


tures such as the optic disc, the blood vessels etc and perfor¬ 
mance of postprocessing steps are well observed from final 
detection result. 

Tabel 2 depicts the result obtained on 89 image each of Di- 
ARETDB1V2 nad EOPHTHAAEX databses. Computation 
of sensitivity and prediction are pixel based. Value of TP, EP, 
TN, EN are measured as number of pixels and eq.(19) is used 
for sensitivity and prediction calculation. Erom the Table 2 it 
can be clearly observed that the presented method is robust to 
other databases also. Eig. 5 shows the variation of sensitiv¬ 
ity vs prediction on the DIARETDB1V2 and EOPhTHAEX 
database. Performance on DIARETDB1V2 with sensitivity 
of 96.54% and prediction of 98.35% is slightly better than on 
EOPHTHAEX with sensitivity of 95.82% and prediction of 
97.85%. Erom this it can be infered that this method is almost 
dataset independent. 

According to Eig. 6, on the DIARETDB1V2 datbase 
AUC value of this method is 0.968 with sensitivty of 96.54%. 
This result is satisfactory for practical purpose of mass 
screening. A comparisons of sensitivity vs 1-specificity 
with recent state-of-the-art method is explored in Eig. 6. The 
method proposed by Zhang et al.[4], achieved AUC value of 
0.95 in the same database. Also the proposed by Giancardo 
et al.[2], achieved AUC value of 0.87. 

Tabel 4 depicts sensitivity, specificity and AUC compar¬ 
isons with other methods. Eventhough not all these meth¬ 
ods have used a common dataset, but this comparions is to 
show the advantages of this method in terms of those mea¬ 
surements. A comparison of AUC values with latest state- 
of-the-art methods is shown in Table 5. Our method perform 
very well in comparison to existing methods. 

Eor classification of soft and hard exudates we have used 
LibSVM Ca package. Each image also converted to Lab 
colour space for getting exudate luminance and colour fea¬ 
tures separately. Training data is a collection of total 5000 
data points having 14 dimension. Training data used both 
DIARETDBV2 and EOPHTHAEX dataset. Eor estimating 
generalisation error a cross validation (10-fold) method was 
used, where training set was divided into 1:9 parts. Eor each 
time 9 parts will be used for training and remaining one part 
for testing. Separate evaluation was done on both dataset for 
soft and hard exudate classification. In Table 3 detailed result 
of classification in pixel level is shown. This final result is 
very satisfactory for practical use of mass screening purpose. 
One of the foremost advantages of our method is that we 
don’t need to detect the optic disc or the vessels separately 



Fig. 6: Sensitivity and 1-specificity variation comparions with 
existing methods on eophthaEX 

to segment out exudates, as this method automatically remove 
those features in detection process. Most of the already exised 
method fH need separate detection for the optic disc and the 
blood vessels in the exudate detection phase. 


Table 4: Comparison results for Exudate detection. Dataset 
are not same 


Methods 

Sen (%) 

Spec(%) 

AUC 

Proposed Method 

96.54 

93.15 

0.968 

Zhang et al. (H 

96.0 

89.0 

0.95 

Niemeijer et al. O 

95.0 

86.0 

0.95 

Sinthanayothin et al. 1161 

88.5 

99.6 

NA 

Osareh et al. [51 

93 

94.1 

NA 

Walter et al. tVfl 

92.74 

100 

NA 


Table 5: Comparison of AUC on two standard database 

Database Proposed 

Zhang et 

Giancardo et 

Method 

al.(2014) 

al(2012) 

eophthaEX 0.952 

0.930 

0.900 

DIARETDBlv2 0.968 

0.950 

0. 930 


4. CONCLUSION 

A novel method for computer aided diagnosis of retinal image 
for exudate detection and analysis is proposed. We have got 
considerable accuracy over two different dataset comprised of 
several varieties of images specifically illuminance changes. 

































with other pathological features etc. Also our machine learn¬ 
ing based classifier SVM works very well for exudate classifi¬ 
cation with designed features. This system can be used for au¬ 
tomated processing of pathological images related to diabetic 
retinopathy, also will be very effective for mass screening. In 
near future, we will incorporate microaneurysms and haem¬ 
orrhages detection to the system to enhance its credibility to 
evaluate the degree of diabetic retinopathy 
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